Engineering correlation and entanglement dynamics in spin systems 
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We show that the correlation and entanglement dynamics of spin systems can be understood in 
terms of propagation of spin waves. This gives a simple, physical explanation of the behaviour 
seen in a number of recent works, in which a localised, low-energy excitation is created and allowed 
to evolve. But it also extends to the scenario of translationally invariant systems in states far 
from equilibrium, which require less local control to prepare. Spin-wave evolution is completely 
determined by the system's dispersion relation, and the latter typically depends on a small number 
of external, physical parameters. Therefore, this new insight into correlation dynamics opens up the 
possibility not only of predicting but also of controlling the propagation velocity and dispersion rate, 
by manipulating these parameters. We demonstrate this analytically in a simple, example system. 

PACS numbers: 03.67.-a, 03.67.Mn 



Correlations play a predominant role in the study of 
spin systems. On the one hand, they characterize dif- 
ferent phases of matter, and thus can help reveal the 
mechanisms underlying phase transitions. On the other, 
they are directly related to the entanglement between 
different spins, which can be exploited by applications in 
the field of quantum information processing. Whereas 
so far, much of the work on correlations has focused on 
the static properties of equilibrium systems, an increas- 
ing interest in the corresponding dynamical properties 
has developed over the last few years. The reason is two 
fold. Firstly, new experimental setups, such as atoms in 
optical lattices, have reached an unprecedented level of 
control, allowing physical parameters to be changed dur- 
ing the experiments. Thus theoretical descriptions of the 
time-dependent properties of such systems have become 
important. Secondly, it has been recognized that the 
way entanglement is created and how it propagates are 
important fundamental questions in quantum informa- 
tion theory. In particular, the answers may influence the 
design of quantum repeaters and networks, whose goal is 
to establish as much entanglement as possible between 
different nodes in the shortest possible time. 

The time evolution of correlation functions in spin 
systems has been studied recently in various scenarios, 
mainly from a condensed matter physics perspective. In 
Refs. 0, 0, two-point correlations were studied numer- 
ically, whereas in Ref. Q their evolution in a critical 
model was studied analytically using conformal field the- 
oretic methods. In all cases, correlations were seen to 
propagate at a finite speed. In Ref. a proof was given 
that correlations necessarily propagate at a finite speed. 
On the other hand, information and entanglement prop- 
agation in spin systems has mostly been studied from a 
quantum information perspective [H, 0, 0, @] • 

In contrast with previous work, we will consider to 
what extent it is possible to control the propagation 
speed and dispersion of the correlations in a transla- 



tionally invariant system, by tuning only simple, global, 
physical parameters. This may be relvant for the op- 
timal creation of entanglement in spin systems, as well 
as contributing to a better understanding of how cor- 
relations are created in dynamical processes, something 
which can be tested experimentally in present setups. We 
will show that, even with this severely limited control 
over the system, the correlation speed can be engineered 
whilst simultaneously keeping dispersion to a minimum, 
so that correlations can be concentrated between partic- 
ular spins. Indeed, by manipulating system parameters 
during the evolution, the speed can be adjusted at will, 
even to the extent of reducing it to zero, allowing corre- 
lations to be frozen at a desired location. 

It is instructive to first consider the entanglement and 
correlation propagation described in the references given 
above from a new perspective. In many of those works, 
correlation propagation can be understood as follows. 
The spin system is initially prepared in its ground state. 
A localised, low-energy excitation is then created (e.g. by 
flipping one spin), and allowed to evolve. Since the low- 
energy excitations take the form of spin waves, the cor- 
relation and entanglement dynamics can be understood 
as nothing other than propagation of spin waves. This is 
completely determined by the dispersion relation (given 
by the system's spectrum). The form of the dispersion 
relation will typically depend on external, physical pa- 
rameters of the system (e.g. the strength of an external 
magnetic field). Thus already in these setups, we can 
manipulate the external parameters to control the dis- 
persion relation, and hence control the propagation of 
correlations. For example, changing the gradient of the 
dispersion relation will change the propagation speed. 

However, the ground state will typically be highly- 
correlated and difficult to prepare, and with the level 
of local control required to create the local excitation 
and break the translational symmetry, more sophisti- 
cated quantum-repeater setups are possible. Also, it is 
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not clear that the correlations will remain localised; they 
are likely to disperse rapidly as they propagate. 

Therefore, we will extend the idea to systems prepared 
in translationally invariant, easily created, uncorrelated 
initial states. For example, the fully polarised state with 
all spins aligned can be prepared by applying a large, 
external magnetic field. As the initial state will be far 
from the ground state, it will contain many excitations. 
The correlation dynamics is then the result of the prop- 
agation and interference of a large number of spin waves 
at many different frequencies. Nonetheless, we will show 
analytically that, at least for some simple models, the 
system can be engineered so that correlations propagate 
in well-defined, localised wave packets, with little dis- 
persion. The external parameters can then be used to 
control the propagation of these correlation packets. 

In the following, we will consider a specific model 
which, despite its simplicity, is sufficiently rich to dis- 
play most of the features we are interested in. The 
model is simple enough to envisage implementing it ex- 
perimentally, for instance using atoms in optical lattices 
or trapped ions. The XY-model for a chain of spin- 
\ particles is described by the Hamiltonian Hxy = 
"jEi((l + 7K^+i + (1 - 7MV+i + 2Aaf), where 
the <t's are the usual Pauli operators and the sum is over 
spin indices. The parameter A can be interpreted as the 
strength of a global, external magnetic field, whereas 7 
controls the anisotropy of the interactions. 

This Hamiltonian can be brought into diagonal form 
by the well-known procedure [9( of applying Jordan- 
Wigner, Fourier and Bogoliubov transformations, giv- 
ing H XY = -|E fe £ fc(7fe7fc ~ 7fc7fc) with spectrum 
e fe = ((cos(2vrfc/7V) - A) 2 + 7 2 sin 2 (2 7 rfc/iV)) 1 / 2 . The 
7^' p are Majorana operators, related to the more usual 
Jordan- Wigner fermionic annihilation operator 7^ by 
7£ = J k + Ik and 7^ = (7^ - Jk)/i, and obey canoni- 
cal anti-commutation relations {7^,7?} = 2Sk,iS a ,b- 

Ultimately, we are interested in "connected" spin- 
spin correlation functions, for example the ZZ correlation 
function Czz{n,m) — (cr^crf„) — (<r*) (of„), in which the 
"classical" part of the correlation is subtracted. These 
are related to the localisable entanglement L(n,m) (the 
maximum average entanglement between two spins n 
and m that can be extracted by local measurements 
on all the others [10|] ) : the natural figure of merit for 
quantum repeaters. In particular, for spin-^ systems, 
L(n,m) > C(n, m) for any connected spin-spin correla- 
tion function C However, we will start by consider- 
ing the simpler, albeit less well-motivated, string corre- 
lation functions such as Sxx{i,j) — (^f(Y\ t< k<j a k) a j) 
(important for revealing "hidden order" in certain mod- 
els [12j|). Their behaviour will give insight into the more 
important spin-spin correlations, and we will use similar 
techniques to calculate both. 

Assume the spin chain is initially in some completely 
separable, uncorrelated state, such as the state with all 



spins down. The interactions are then switched on and, 
as this initial state is not an eigenstate of the Hamilto- 
nian (unless A — > 00), the state evolves in time. The 
initial state is the vacuum of the Majorana operators 
x i{pi) = Ylj^j^i obtained after applying just the 
Jordan- Wigner transformation, and is completely deter- 
mined by its two-point correlation functions. In other 
words, the vacuum is a fermionic Gaussian state, and 
can be represented by its covariance matrix r m „ — 
\ ([r m ,r n ]) where r 2 /_i = xi and r 2 i = p/. 

From the Heisenberg evolution equations, it is sim- 
ple to show that any evolution governed by a quadratic 
Hamiltonian corresponds to an orthogonal transforma- 
tion of the covariance matrix. It is also clear that, as 
the Fourier and Bogoliubov transformations are canon- 
ical (anti-commutation-relation-preserving) transforma- 
tions of the Majorana operators, they similarly leave 
Gaussian states Gaussian, and they too can be expressed 
as orthogonal transformations. Thus the time-evolved 
state of the system is given by a series of orthogonal 
transformations of the fermionic vacuum: 

r(t) = or vac o T , o = o£ T o£ og o(t). (i) 

This is a block- Toeplitz matrix, composed of 2 x 2 blocks 
G x at distance x from the main diagonal. In the thermo- 
dynamic limit N — > 00 with — > <fi and Ek — > s(4>) = e, 

G x = [ d(f> ( 90 91 ] , 5o = iSsm(<f>x) sin(2ei) 
J-t; \9-i 9oJ 

g±x = 2CSsin(<j)x) sin 2 (rf) ± cos{4>x)(c 2 + S 2 cos(2rf)) 

where C = (cos(</>) — X)/e((f>), S = 7 sin(</>)/£(</>), and 
x = m — n. We can now calculate certain string corre- 
lations, which are given directly by elements of the co- 
variance matrix. For example, (^(n„< l<m O^) = 

7r 2 „_l,2m-l = Z) a =±l s I-* d^ 5 " COS((/>X + 2set)/2. 

Although the evolution of the string correlations is pro- 
duced by the collective dynamics of a large number of 
excitations, this expression has a simple, physical inter- 
pretation: it is the equation for two wave packets with 
envelope S/2 propagating in opposite directions along 
the chain, according to a dispersion relation given by the 
system's spectrum e((j>). This wave-packet interpretation 
allows us to make quantitative predictions as to how the 
dynamics will be affected if the system parameters 7 and 
A are modified. Specifically, modifying the parameters 
will change the dispersion relation, changing the group 
velocity of the correlation packets, as well as the rate at 
which they disperse. (The wave-packet envelopes also de- 
pend on the system parameters, so the relevant region of 
the dispersion relation may also change.) Thus by vary- 
ing only global physical parameters, we can control the 
speed at which correlations propagate. 

Does this hold true for the more interesting spin-spin 
correlations? We will show analytically that they have 




a similar wave-packet description, although in terms of 
multiple packets propagating simultaneously. This will 
allow us to predict the behaviour of the spin-spin corre- 
lation dynamics for different values of the system param- 
eters. In particular, we will show that the correlations 
can be made to propagate in well-defined packets whose 
speed can be engineered by tuning the system parame- 
ters. Moreover, the propagation speed can be controlled 
as the system is evolving, so we can speed up or slow 
down the packets, even to the extent of reducing the 
speed to zero. We confirm our predictions by numeri- 
cally evaluating the analytic expressions. 

Let us now calculate the spin-spin connected correla- 
tion function using the covariance matrix derived above. 
We have ct* = x n p n , so the ZZ connected correla- 
tion function is given by C Z z{x) = (x n p m ) (p n x m ) - 
(x n x m ) (pnPm), where we have used Wick's theorem to 
expand the expectation value of the product of four Ma- 
jor ana operators into a sum of expectation values of pairs 
[I3|,ll4(. The latter are given by covariance matrix ele- 
ments, resulting in the following analytic expression for 
the correlation function: 



C' zz (x,t) 2 



(2) 
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(scos(0x — 2sef)) 

2 s=±l 

d())CS(sm((j)x) sm(<px + 2set)) 

2 s =±i 

dej) (C 2 cos(cj)x) + — ^2 c °s(</>£ + 2set)) 
2 „=±i 



Although more complicated than the string correla- 
tions, this expression also describes wave packets evolv- 
ing according to the same dispersion relation e(<fi), al- 
beit multiple packets with different envelopes propagat- 
ing and interfering simultaneously (three in each direc- 
tion). In many parameter regimes, broad (in frequency- 
space) wave packets and a highly non-linear dispersion 
relation will cause the correlations to rapidly disperse 
and disappear. However, we can find regimes in which 
the wave packets are located in nearly linear regions of 
the dispersion relation, and maintain their shape as they 
propagate. For example, at 7 = 1.1 and A = 2, all three 
wave packets of Eq. j2]) are nearly identical, and reside 
in an almost-linear region of the dispersion relation with 
gradient roughly equal to 2, as shown in Fig. [T] (inset). 
The spin-spin correlation dynamics will therefore involve 
well-defined correlation packets propagating at a speed 
dx/dt — 2, dispersing only slowly as they propagate. 
Fig. [2 shows the result of numerically evaluating Eq. @, 
which clearly confirms the predictions. 

We can engineer a different correlation speed by chang- 
ing the parameters. For instance, for 7 = 10 and A = 0.9 
we predict a higher propagation speed dx/dt ~ 18, al- 




FIG. 1: For 7 = 1.1, A = 2, all the wave-packet envelopes 
from Eq. ([2]) (non-red curves, inset) are similar in form, cen- 
tred around a nearly linear region of the dispersion relation 
with gradient ~ 2 (red curve, inset). Thus the correlations 
C zz (x,t) (indicated by the shading, main plot) propagate in 
well-defined packets at a speed given by the gradient. 



though at the expense of increased dispersion. The nu- 
merical results of Fig. [5] show precisely this behaviour. 
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FIG. 2: For 7 = 10, A = 0.9, the wave-packet envelopes 
are spread over the entire frequency range (non-red curves, 
inset). However, the dispersion relation (red curve, inset) is 
almost linear for wavenumbers not too near 7r, with gradients 
±18. As the envelopes are symmetric about n, most of the 
correlations C zz {x,t) (shading, main plot) still propagate at 
a well-defined speed « 18, faster than in Fig. \T\ (note scale), 
though as expected they also show more dispersion. 

An even more interesting possibility is controlling 
the correlation packets as they propagate. If the sys- 
tem parameters are changed continuously in time, the 
XY-Hamiltonian becomes time-dependent, and the or- 
thogonal evolution operator 0(t) in Eq. (fT]) is given 
by a time-ordered exponential 0(t) = T[e^> dt )] = 
lim/i_o ni=i"' e A ( nh \ (A is a time-dependent, anti- 
symmetric matrix determined by the Hamiltonian.) In 
general, the time-ordering is essential. But if the system 
parameters change slowly in time, dropping it will give 
a good approximation to the evolution operator. The 
state at time t is then just given by evolution under the 
time-average (up to t) of the Hamiltonian. If we remain 
in a parameter regime for which the relevant region of 
the dispersion relation is nearly linear, adjusting the pa- 
rameters changes the gradient without significantly af- 
fecting its curvature or the form of the wave packets. 
Thus, to good approximation, slowly adjusting the pa- 
rameters should control the speed of the wave packets as 
they propagate, allowing us to speed them up and slow 
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them down. Numerically evaluating the time-ordered ex- 
ponential shows this is indeed possible (Fig. [3]). 
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FIG. 3: Starting from 7 = 1.1, A = 2 as in Fig. Q] 7 and 
A are smoothly changed to move from the red to the green 
dispersion relation (inset), increasing the correlation speed. 

Clearly it would be useful to be able to stop the correla- 
tions once they reach a desired location. One way would 
be to simply switch off the interactions. But strictly 
speaking this would require more control than is provided 
by the two parameters defined in the Hamiltonian (there 
is no value of 7 for which all interaction terms vanish), 
and may be difficult in physical implementations. If the 
spin model were realised in a solid-state system, for ex- 
ample, switching off the interactions would likely involve 
fabricating an entirely new system. In any case, we will 
show that switching off the interactions is not necessary 
in order to freeze correlations at a specific location. 

Instead of changing the parameters continuously, we 
now consider changing them abruptly. The time-evolved 
covariance matrix in this scenario can be calculated an- 
alytically by the same methods as used above. Suppose 
the initial system parameters 70 and Ao are suddenly 
changed to 71 and Ai at time t\. The spin-spin correla- 
tions will initially evolve according to Eq. 0, as before. 
After time t\, the evolution becomes more complicated. 
The analogue of Eq. ^ separates into a sum of wave 
packets evolving in four different ways: those that ini- 
tially evolve according to Eq and subsequently (after ti) 
evolve according to e±, those that subsequently evolve 
according to — e\, those that only start evolving at t\, 
and those that that undergo no further evolution after 
ti (Fig. |4} . For t > tt, the terms whose evolution is 
"frozen" at time t\ are given by 




Since t does not appear on the right hand side, this ex- 
pression clearly describes wave packets that propagate 
until time t\ and then stop. Using these, we can move 
correlations to the desired location, then "quench" the 



system by abruptly changing the parameters, freezing the 
correlations at that location, as shown in Fig. [4] 
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FIG. 4: The system is initially allowed to evolve with 70 = 0.9 
Ao = 0.5, then "quenched" at time ii = 20 to 71 = 0.1, 
Ai = 10. Some of the correlations are frozen at the separation 
(x w 20) they reached at £1. Others propagate according to 
the new dispersion relation, or are "reflected". 

We have shown that entanglement and correlation 
propagation in many spin-model setups can be under- 
stood in terms of propagation of spin waves, and have 
introduced the idea of controlling the dynamics via their 
dispersion relation, by manipulating external parameters 
of the system. Although this is in principle possible 
for almost any spin system, preparing a single-excitation 
initial state would require control over individual spins. 
Therefore, we have analysed in detail the more complex 
case of systems prepared in uncorrelated, translationally 
invariant initial states, which typically contain many ex- 
citations. We have shown for an example model that the 
dynamics can be described by a small number of corre- 
lation wave packets, and that the control afforded by a 
few external, physical parameters is sufficient to allow 
detailed control over the propagation of correlations. 
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